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Let a,b,c be positive integers with b< min (a,c— 1). In an infinite sequence of indepen- 
dent trials with different success probabilities, define a "run" to consist of a consecutive 
successes, and a "quota" to consist of c or fewer consecutive trials containing b or more 
successes. Recursion formulas are derived for the probabilities governing the first appear- 
ance of a run immediately preceded by a quota. 

1. Introduction 

This paper deals with the following problem, which arises in the analysis of certain target- 
detection systems. Consider an infinite sequence {T n : n=l, 2, . . . } of independent probabilistic 
"trials." T n yields a "success" with probability p n , and a "failure" with probability q n =l —p n > 
A run is defined to be a set of a consecutive trials all yielding successes, where a is a positive 
integer specified in advance. A quota is defined to be a set of c consecutive trials, of which al 
least b yield successes; * here b and c are positive integers specified in advance, with 6<V. 
When does the first run which immediately follows a quota occur? 

Of course the question can only be answered on a probabilistic basis, i.e., one seeks the 
probability, as a function of n, that T n is the last trial of the first such run. The analysis 
given below leads to a recursion equation for these probabilities (more precisely, for certain 
other closely relate:! ones) under the additional hypothesis that a>b. The equation is some- 
what unpleasant looking, but offers no particular difficulty as regards implementation on a 
digital computer. 2 Both the analysis and its result appear likely to be more complicated if 
a<Cb, so that it is fortunate that the hypothesis a>b covers the applications we have in mind. 

The special case a = 2, b = l, c=2 has been solved by E. Brookner as part of a larger study ; 3 
the first author appreciates access to an advance copy of Brookner's paper. Since the methods 
used below are rather ad hoc, it should be noted that other work by Brookner 4 provides, in 
principle, a systematic approach to such problems; in particular, his methods can apparently 
be used to derive a solution in the case a<Cb. 

The next section contains a precise statement of the problem, and a guide to the main 
results of the investigation. 

2. Notation and Summary of Results 

The probabilistic events of interest will be regarded as subsets of a sample space, 5 and so 
the notation of Boolean set-algebra will be employed. Thus the symbols U and f] stand for 
set-theoretic union and intersection, respectively, while A' will generically denote the com- 
plement of the event A. We depart from standard notation by using the symbols S and "-)-" 
for the set-theoretic union of disjoint events (i.e., "events whose probabilities can simply be 
summed"). We set p n =0 for n<^l. 



*Sponsorcd in part by the North American Air Defense Command (Directorate of Operations Analysis) and U.S. Army Air Defense 
Command. 

1 This definition must be slightly modified near the beginning; of the sequence; see the definition of Q n in section 2. 

2 This assumes that a routine for computing the sequence \p n \ is available. 

3 Eli Brookner, Cumulative probability of target detection for pulsed surveillance radars, Columbia University Electronics Research 
Laboratories Tech. Rept. T-10/122. 

4 Eli Brookner, Determination of the statistics of a recurrent event in a Markov chain, Columbia University Electronics Research Laboratory 
unpublished manuscript. 

5 For this approach to probabilistic events, see e.g., the'standard reference by W. Feller, An introduction to probability theory and its ap- 
plications (John Wiley & Sons, Inc., New York, N.Y.). 
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A set-theoretic union over an empty set of indices will be taken, by convention, to be the 
null event 4>; similarly a numerical sum over an empty set of indices will be taken as zero. 
A set-theoretic intersection over an empty set of indices will be taken as the "sure event" 
(i.e., the entire sample space), while a product over an empty set of indices will be assigned 
the value unity. 

The following events will be involved in the analysis : 

&n'- T n yi elded a success {S n =--(j> for n<^l), 
S' n : T n yielded a failure, 

a-l 

R n = fl S n -j'. T n is the last trial of a run (of a consecutive successes 6 ) 
o 

For the exact definition of "quota," it is convenient to define 

m(w)=max (n— c+1,1). (1) 

Then we have the definitions 

Q n : at least b trials among T n , T n _i, . . ., T m(n) yield successes, 
E n =B n [}Q n - a : the type of event whose first occurrence interests us. 

Verbally, E n occurs if the a trials numbered n, ?i— 1, . . ., n— a+1 all Afield successes, and in 
addition b or more of the min (c,n — a) trials numbered n — a, n — a—1, . . ., m(n — a) yield 
successes. Still other events will be required and introduced later. 

The generic notation P(A) will be used for the probability of an event A. The probabili- 
ties of interest, those referring to the first trial which ends a run immediately preceded by a 
quota, are therefore designated 

z,=p(^nV^;.) (2) 

It will be simpler to deal, instead, with the probability 

r,=p(n #;) (3) 

that no run preceded by a quota occurs up to and including the trial T n . The sequences 
{X n } and {Y n } are related by 

n 

Y n —\ — S X i} X n =Y n _ 1 —Y n , (4) 

i=i 

and the sequence {Y n } clearly obeys the "initial conditions" 

Y n =\ \in<:a+b, (5) 

a+b 

Y a+i =i- n Pt . (6) 

i = \ 

We therefore assume n^>a-\-b in all that follows. 

The reader is exhorted to aid himself, in following the derivations, with the diagram of an 
axis on which points corresponding to T ri} T n _ a+1 ,T n _ a , T„_ a _ c+W , T n _ a _ c+b and T n _ a _ c+l are 
marked off. 

The ultimate goal is a recursion equation for the sequence {Y n }, while the method is the 

11 

decomposition of fl E t into a union of disjoint events, each the intersection of certain inde- 
pendent events. 7 This analysis is carried out in section 3. The desired recursion equation, 



e This run may be imbedded in a still longer sequence of consecutive successes. 

" If the set of trials whose outcomes determine whether or not an event A occurs is disjoint from the corresponding set of trials for event B, 
then A and B are clearly independent. This sufficient condition for independence is the one implicitly appealed to throughout the paper. 
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for a-\-b<^n<2a+c, is derived without (he assumption a>b and is given in eq (24). For 
ny>2a+c, under the assumption a>b, it is given in eq (27). Both of these equations involve 
as coefficients certain auxiliary probabilities whose calculation is the subject of section 6. As 
for the case n^>2a-\-c and a<^b, which is not solved in this paper, our best progress to date is 
reported in eq (31). 

The assumption a>b is made throughout sections 4, 5, and 7. In section 4 it is shown 
that the original recursion equation can be replaced by one which has higher order but fewer 
nonzero terms, and so may be more suitable for computation. The new equation is given 
by eq (34), which is valid for n>2a + c and 6>1, by eq (34a), which is valid for a-\-b<Cn<C2a + c 
and 6>1, by eq (34b), which is valid for n^>2a+c and b = (), and by eq (34c), which is valid for 
n^>2a J rc and 6=1. The remaining cases (n = 2a J rc, or n<^2a-\-c and 6<1) are left to the 
reader. 

Section 5 contains an examination of the special case in which p n is independent of n. 
The general recursion eq (24) becomes eqs (37) and (38), while the general recursion eq (27) 
becomes eq (39); for 6>1, the general eq (34) yields eq (40). For this special case (i.e., p n 
=p and 6^>1) a method is given for calculating the mean occurrence time for the first run 
preceded by a quota; this is given compactly as eq (47). The extension to b<\ offers no 
difficulty. 

In section 7 we study the probability Z n thai the first run preceded by a quota, if it ends 
at T nj occurred as a sequence of a-\-b consecutive successes. (Note that this is not the same 
as the probability that the first sequence of a + 6 consecutive successes ends at T n .) Equation 
(57), which is valid for //>a + c+l and 6>0, shows how the sequence {Z,,} is determined by 
the sequences \X^ and {F n }. The illustrative case 6=1, c = 2 is investigated in some detail, 
and bounds on Z n are found in terms of p n ; see e.g., eq (61), valid for ?/>3a-f 3 and a>l, 
and eqs (62) and (63), valid for //>(> and a=l. 

3. Analysis 

In this section the hypothesis a>b is not at first assumed. From E n =R n f] Q n -a, one has 

E' n =R' n +(R n nQ' n - a ), 

where the convention for "disjoint union' 7 has been used. Thus 

n E'i=E' n {\ V #;=(#; n VsjV(fi,n&-.n n n e'\ (7) 

i=l i=l \ i=l / \ t=l / 

The first disjunct on the right in eq (7) will now be decomposed further. Let R nk be the 
event that the first failure among trials T n ,T n ^ 1) . . ., (in that order) occurs at T n - k . Then 

a-l 

#»=2 R'nic, which implies that 8 



A=0 



n — 1 a — I / n—1 \ 

Kn n £5=s(fli*n n e'X (8) 

i = l A-=0 \ i=l / 

However, it is easy to show the inclusion 

R'nk^E'i for n-Jc<i<n-l, l<k<a-l; (9) 

for if R nk occurs then so does S n - k , and obviously this precludes the occurrence of R t (and 
thus of E f ) for n — k<i and k<a— 1. By (8) and (9), we have the desired decomposition 8a 



s The distributive law of Boolean set-algebra is to be used. It will be employed at other points without explicit mention. 
8a Note the use, for k=0, of our convention about intersections over an empty set of indices. 
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n — 1 a— 1 / w— 1 »— fc— 1 \ 

tf;n n sKS(i?i*n n ^ n s 

i = l fe=0 \ i=ra-fe »=1 / 

a — / n— k— l \ 

=S(fli*n n #;)• ( 10 > 

fc=0 \ i=l / 



7i— fc-1 
Now i?^fc and fl 2££ are independent, and furthermore 

i=l 



fe-1 /n-k-1 \ 

p(iC)=<z»-*n #»_.,, p( n bji-f^-i. (ii) 

;=0 \ 1=1 / 



fe-1 /n-k-1 

J 

Equations (7), (10), and (11) therefore yield 



a-l /k-1 \ / n-1 \ 

r B =s&-» n^ 7,.^+p B n n<2Un n e'A (12) 

fc=0 \j=0 / \ i = l / 

Equation (12) indicates that further progress requires a decomposition of the event 

n-l , 

Q' n .. a f) fl E t , As a first step in this direction, it will be shown that 

2 = 1 

Qn-a^Ei foTn-a<i<n-l, (13) 

which implies that 

n—l n—1 n—a—l n—a—l 

Q'„-aC\ n #;=<?;_« n n #;n n #;=#;-„ n n e\ 

i = l i = n—a i = l i = l 

and thus that 

Rn fl Q'n-a fl ? 'fl' #{ = #» fl Q'n-a fl " H ' S{. (14) 

t=l »=1 

, , a— 1 

Since (a) R n is independent of (? n _ a and also of E t for i<n—a, and (b) P(R n )— H p n _j, 

j=o 

application of eq (14) to eq (12) yields 

o-l /k-l \ /a-l \ / n-a-1 \ 

Y n =J2 a-* ( n Pn -j ) Yn-t-M n Pn _j )P{ Q'n-a [\ n s; )■ ( 15 > 

£ = \j=0 / V=0 / \ t=l / 

To prove (13), assume that Q n - a and E % both hold with n—a<i<n — l. Then R t holds, so 
that i>a+6 and the trials T i} Ti^ h . . ., TVa+i all yield successes; in particular, this holds for 
the trials 

J-n-ay J- n-l-aj • • •> -* i+l-a* (16) 

If %<n— b, then T n - a ,T n - a -hT n -a-2, • • ., T n -. a -b+i would be in the sequence (16) and there- 
fore would all yield successes, contradicting the occurrence of Q' n - a . So ^>7l— 6, i.e., 
i—a^>n—a—b^>n — a—c, a fact to be used below in (17), in the form i—a>m(?i—a). 

Since E t occurs, Q t ^ a must occur, so that there are at least b successes in the trials 
Ti-arTi-a-u • • •> T m{ ^ a) . This last-mentioned sequence of trials can be split into two 
subsequences, 

-*- i— a) J- i—a— 1? • • •> J- m(n—a)j \1 '/ 

and 

■*■ n— a— o -L n— a— c— lj • • •? -*- m(i—a) v-*-*/ 

where the second subsequence is absent if n < a+c. Of the "at least b successes" just mentioned, 
at most n—i can occur in the subsequence (18), since the number of trials in this subsequence is 
at most 

n— a— c— (i— a— c+ 1) + 1 =n-~ i. 
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Thus there are at least b—(n—i) successes in the subsequence (17). Combining these with 
the n— i successes in the sequence (16) (which is disjoint from (17)) yields at least b successes 
in the sequence T n - a ,T n _ a _ 1} . . ., T m(n - a) , but this contradicts the occurrence of Q n -. a > Thus 
the assumption that Q n - a and E t both hold is untenable, and the inclusion (13) is proved. 

Equation (15) shows that further progress requires a decomposition of the event 

/i— a- 1 , 

Q' n - a D n E { . This will first be done for n<2a-\-c i by showing that 

i = \ 

Q a -nEL E i for l <i<n—a — l, if a+b<n<a+c, (19) 

Qn-aEE'i f or n— a— c+l<i<n— a— 1, if a+c<n<2a+c. (20) 

These results imply the decompositions 

Q' n - a (] fl K'i = Q' n - a if a+b<n<a+c } (21) 

i = l 

r^_ a n" n sj^^.n" n e E' t a a+c<n<2a+c. (22) 

Since £,\ cannot occur for ?<^, and can occur for i=a only when 6 = 0, in which case we adopt 
the convention c=0, the last two displays can be rewritten as 

Qn-ad 71 n l ^; = ^;_ a if a+6<w<2a+c. 

Substitution of these results into eq (15) yields 

r w =§ &-* f n />„_,) Y^^+Cn p n _i) P(Q' n . a ) for a+6<><2a+c. (24) 

The evaluation of P(Q' n _ a ), which appears in eq (24), will be discussed later. 

To prove the inclusions (19) and (20), first note that both can be simultaneously written as 

Qn-cS^E't for m(n— a)<i<n— a— 1, if a+b<ji<2a+c. 

Suppose, on the contrary, that both Q n - a and E t occur. Since E t occurs, there must be at 
least a-\- b successes in the sequence of trials T t ,Ti- lt . . ., T m{i _ n) . Since 

i<n— a— 1 <(2a+c)— a— l=a+c— 1 

we have m(i— a) = l, and so there are at least &+& successes among the trials 7 7 i ,7V I , . . ., TV. 
Since Q w _ a holds, at most 6—1 of these successes occur among trials Ti,Ti- h . . ., T m0i _ a) , 
and so there are at least a+b — (b — l)=a+l successes in trials T w(n _ a) _ 1 , T mU _ a) _ ? , . . ., 7\. 
Thus 

a+1 <m(n— a)— l=max (n— a— c,0), 

implying a+l<n—a—c and thus n^>2a+c J a contradiction. So (19) and (20) are proved. 
Thus the desired recursion equation for the sequence {F n } has been derived for n<2a+c, 
and is given by eqs (5), (6), and (24). From now on, n^>2a J rC is assumed. 

it— a — 1 

To decompose Q' n - a r\ fl E\ when n>2a+c, it is helpful to observe that Q' n - a can be de- 

i=\ 

fined as the occurrence of at least c— (b — l) = (c — 6) + l failures in the c trials between 7 7 „_ a 
and T n _ a _ c+1 inclusive. Not all of these failures can occur in the c—b trials between T n _ a 
and T w _ a _ (c _ 6)+1 inclusive. We therefore introduce the new events defined for 0<£<6— 1 l>\ 

F„,: the (c — 6) + 1st failure among T n _ a , T n _ a _ h . . ., 

Tn-a-c+i (in that order) occurs at I 7 w _ a _ (c ._ 6) _ r . 

6-1 

Then Qn- a =^2 F nt , which implies the decomposition 
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n—a—1 6—1 / 7i— a— 1 \ 

<?Un n i?;=s(F B( n n e'A, 

i=l t=0 \ i=l / 

so that eq (15) becomes 

fl-i /k— l \ /«-i \ &-i / w-fl-i \ 

F„=s &_* ( n />„_,. ) !%,_,_,+( n /;„_ ; ) s n f»i n n #; )• (25) 

The solution under the extra hypothesis a>b will now be completed by showing that 

FntEiE't for ^-a-(c-6)-/<i<//-a-l ; if a>b. (26) 

This will imply 

n — a — 1 /? — ft — (c— 6) — £ — 1 

! = 1 1=1 

which, since 2^, and E' t are independent for l<i<n—a—(e—b)—t—l, will turn eq (25) into 

r»=2 a,_» ( n 1 p..,) TVt-.+fn p,. i )sP(f„)r M .„. w . ( . 1 if «>6. (27) 

Equation (27) is the desired recursion relation for n^>2a-\-c; together with eq (24) and the 
initial conditions of eqs (5) and (6), it entirely specifies the sequence (F w ) when a>b. The 
calculation of the coefficients P(F nt ) will be discussed later. 

To prove the inclusion (26), suppose that F nt and E t both hold, with a>b and 
n—a—(c—b)—t<i<n—a—l. Since i>n—a—(c — b) — t, the (success-yielding) trials be- 
tween T t and T ? :_ a+ i inclusive must have greater subscripts than the (failure-yielding) 
T n - a -(c-b)-t, so that 

n— a— c+1 <n—a— (c—b)—t<^i—a-\- 1 < i<n— a— 1 

so that the a successes just mentioned fall between T„_ a and T w _ a _ c _ f _ l inclusive. Because 
a>b, this implies that Q n - a occurs, contradicting the assumption that F ni occurs (recall that 

Fn&Qn-a). 

The reasoning just employed carries over in part to the case a<b. Specifically, it will be 
shown that 

F nt ^Ei for n-a-(c-b)-t<i<n-a-l and*<a-l (28) 

so that eq (25) becomes 

/k-l 



a-l /k-l \ /a-l \ a-l 

F,=Sfr-4 nPn-j)Y n _ k M n Vn-AHP{Fnt)Y n - a - 
fc=0 \;=0 / \;=0 / t=0 

(a-l \ b-i / n—a-1 

nft-, S? W n /<:;). (29) 
;=0 / *=a \ i=l / 



-(c-b)-t-l 

h-l 



To prove the inclusion (28), suppose that F nt and E t both hold. As in the last paragraph, 
E t requires at least a successes between trials T n _ a and T n _ a _ ic _ b) „ t inclusive, whereas F nt 
requires exactly (c— 6) + l failures (and thus exactly t successes) in these (c — 6)+^+l trials. 
Therefore, a<t, contrary to hypothesis in (28). 
It is also easy to show that, for a<t<b— 1, 

F nt SE'i iovn-a-(c-b)-t<i<n-(c-b)-t-l 7 (30) 

so that eq (29) becomes 



a-l /k-l \ /a-l \ a-l 

k=0 \;=0 / \j=0 / *=0 

(a — 1 \ 6—1 / n—a—1 n-a-(c—b)—t-l \ 

n?,- ; S'' n,n n #;n n £;). (3i) 
j=0 / t = a \ i = n-(c-b)-t i = l / 



-a- (c-6) -1 — 1 
6—1 / n—a—1 n—a — (c—b)—t — l 
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for if F nt and E t bolh occur, with \>n— a— (c— b)— t, then the fact that T u 7 7 ,-_n . . ., 2Va+i 
all yield successes whereas T n - a -{c-T>)-t yields a failure, implies that 

i—a-\-iy>n—a—(c—b)—t, 

i.e., i>n — (c — b) — t. The authors have not succeeded in carrying this analysis further in 
the case a<Cb. 

4. A Possible Simplification 

The recursion ea (27) is of order a+c, i.e., each member of the sequence {Y n } is expressed 
in terms of the a+c members immediately preceding it. Actually, only a-\-b of these a + c 
preceding members are employed, 9 the others appearing with zero coefficients in eq (27). In 
this section eq (27) is used to derive a recursion equation of the higher degree a + c+1, which 
however has the advantage that calculation of a particular member of {Y v } requires only 
6 + 3 of the preceding a + c+1 members (and only 2 if 6=0). Of course a detailed evaluation 
of this "advantage" would involve a comparison of the computational labor for the coefficients 
of eq (27) versus that for the coefficients of the new equation. 

For the derivation, first change n to // -f- 1 in eq (27) and shift "dummy indices" to write 
t he result , for 6 1 , as 

Y*+l = qn+lY H +^ q.n-k{ H p n +l-j ) F.-jt-l + f ^ Pn+l-j)( U q n + l-a-i) Y n-a-(c-!» 

+(nV+w)P(ft+i.i)^ (32) 

Next multiply eq (27) through by p n +i, < l,1( l write the result as 

0=p n + l Y n — ]TJ 2»-*( n Pn+l-j ) Yn-jc-t— q n -a+l ( TL Pn + l-j) Y 'n-a 

— ( II Pn+l-A ( II Zn-a-i) Y n-a- (c-») - l — ( N Pn+l-A ZJ P(F Ht )Y H -a- (*-&)-*- 

Now add eqs (32) and (33), noting that £> w+1 +</ n+1 =l: 

^ « 4 i = ^ »— Q.n-a+1 ( n p n +i-j) I «- a + ( H p n +l-j ) ( H </ w + i-„-/ j } „_ a _ Cc _ &) 

+(n p n + w J P(^»+i f i)— ^*+i-« (n g«- a _ J r„ 
+(nV + i-y) 2 [p^+li+O-a+i-^^i)] F«. 

-fnA+l-/) P^^-OI^-a-c (34) 

This is the desired equation. 

A similar equation can be derived to replace eq (24). The procedure is as above, and the 
result is: 

Y n ^ = Y n -g n _ a+l (Vpn+i-y) Y n - a +(*j£p n+l -?\ [P(Ci+i-«)"-A+i-«P(Ci- a )] (34a) 
fora-f-6<//<2a + c. The analogous equation for the "transition case" n=2a-\-c will be omitted. 

9 The order of the recursion is related to the "memory'' requirements of a possible computer program, while the number of nonzero coefficients 
is relevant to the amount of computation required for each Y n . 
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(33) 



-a-(c-b)-l, 



-a-{c-b)-t-l 



If 6 = or 6=1, the preceding equations must be modified. For 6 = 0, eq (34) becomes 

Y n +i=Y n — q n - a +\ ( n Pn+i-A Y n -a, (34b) 



while for 6=1 it becomes 

/c-b 



Y n+ i = Y n —q n - a+ i I II Pn+l-j ) Y n - a +l II jp»+i_j( II g n+1 _ a _ f j !F n . a _ c+1 

-(jft+J ( H Zn-a-i} Yn-a-o (34c) 

The analogous results for eq (34a) are omitted. 

5. A Special Case 

It is natural to inquire what form the recursion equations take in the relatively simple 
case when all the trials are governed by the same probability distribution, i.e., 

Pn=P, 2»=2=1— P- 
The initial conditions are 

Y„ = l for n<a+b, (35) 

Y a+h =l-p a + b . (36) 

From the definition of Q n ~ a it follows that in eq (24) 

P(Q:-a)^( n ~ a: )p t (l n - a - t for a+b<n<a+c, 

P«?U)=S(5)^2 c " r for a +c<n<2a+c. 

Thus the early recursion equations read (for q^>0) 

Y n =q sV-i-i+rt"" 1 S f^V^"' if a+60<«+ c , (37) 

I>gg^V.-i+rt c (g (S)^" 1 ) Yn-a-c if a+c<»<2a+c. (38) 

From the definition of F nt , it follows that 

P(^«)==( c ~* +f )pY- 6+1 , 

so that eq (27), which is valid if a >b, becomes 

?»=? S/r„_ t .,+| ) Y-»+' E ('"J+^T,.,.,.,.,., iin>2a+c, (39) 

while eq (34) becomes (for 6>1) 

r„ +1 =r re -^r 7! _ a +pY _6+1 ^- a -( C - 6 )+(c-6)^ +1 2 c - 6+1 F B _ a _ (c _ w _ 1 

+i> a +v -» + i g ^-|+ < )yr„_ a _ (c _ S) _ / _ 1 - 2 )«+v-»+' (J-|) r n _ a _ c . (40) 

Let the random variable A 7 be defined by the condition that T N is the last trial of the 
first run immediately preceded by a quota. The sequence {X n } defined in eq (2) is thus 
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precisely the frequency distribution function of TV. Although the equations above are not 
sped acularly simple despite the simple choice of the sequence {p n }, they at least permit descrip- 
tion of a finite procedure for calculating the expected value E(JST) of TV in the special case under 
consideration (with a > b). 

To this end, first use eq (4) to write 

E(N) = S nX n =X, +X>(F„ _ ,-¥„) 

1 2 

£ nY n ^-JZ nY n ) 

2 2 / 

=l+lim (T,Y n -mY,\ 

Since the definition (3) of the sequence { Y n ) ensures that this is a monotone nonascending 
sequence of noimegative terms, a classical theorem on infinite series 10 can be invoked to assert 

that mY m -*0 if 2 F n <> . This yields 

i 

E(N)=l+T,Y n ifSI r »<«, (41) 

i i 

so that the evaluation of E(N) reduces (o that of ]C Y n . This of course holds even in the general 

i 
case. 

We now employ the standard formal procedure for summing an infinite series whose terms 
obey a linear difference equation with constant coefficients. Write eq (39) briefly as 

F,=S atYn-t for n>2a+c, (42) 

i = \ 

and introduce the generating function 

5(2)=SF,2" (43) 

1 

as well as the polynomial 

/(*) = 1-!8<M'. (44) 



i=l 



An elementary calculation using eq (42) yields 



1 2 = 1 1 

1 i=l « = i-f-l 



=S^s n -2 (s n S {^F.-! : l<i<min(a+c,n-l)}) 

1 n=2 

a+c / n-l \ / a+c \ 

= 1+S ( F.-S «<F„_ 4 ) S"+ S ( F.-S "iYn-t ) :" 

2 \ t=l / a+c+1 \ i=l / 

a+c / w-1 \ 2a+c / a+c \ 

=1+S (F.-S a,F,_,) z»+ 2 (r,-SaiF._,)a»=A(8), (45) 

2 \ 2=1 / a+c+l \ / = ! / 



where the coefficients of A(s) can be found from eq (35) to (38). 
It will be shown that for 0<p<l 



f (2)9*0 forjs|<l. (46) 



"> See p. 124 of Knopp, Theory and Application of Infinite Series, Blackie (1951 edition). 
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This justifies the preceding formal manipulations for |2|<1, permits solution of eq (45) as 

g(z) = h{z)lj(z), and permits use of eq (43) to write ^,Y n =g(l), so that eq (41) becomes 

i 

E(N)=l+h(l)/f(l), (47) 

the promised formula for E(A r ). In connection with eq (47), note that by eqs (39) and (44) 

f(\) = p a (l-cr' +1 'S( C ~ b t + t )p') (48) 

To prove the nonequality (46), note that all a t in eqs (42) and (39) are nonnegative, 
so that for |s|<l, 



a-\-c 






a-\-c 

thus (46) is true if XI ^<1, i.e> ? if /(1)>0. By eq (48) this is equivalent to 



§n + >¥-«<i, 



which is true since the left-hand sum is the value at b — 1 of the cumulative negative binomial 

probability distribution with parameters c— 6 + 1 and q. (See Feller, oj). cit. in footnote 5.) 

As in section 4, the cases 6=0 and 6=1 require special treatment, which is left to the reader. 

6. Auxiliary Probabilities 

This section deals with the probabilities P(Ql_ a ), which appear in eq (24), and P(F nt ), 
which appear in eq (27). The "best" way of calculating these quantities (or more generally, 
of arranging the whole recursive computation) is likely to depend strongly on the values of 
a,b,c } on the peculiarities of the specific sequence [p n ] of success probabilities under considera- 
tion, and perhaps on special characteristics of the calculating equipment used. Thus the pro- 
cedure sketched below is illustrative, rather than " recommended." 

For a+6<?i<a+c, Q' n -a * s the event that the sequence of n — a trials T n _ a ,T n _ a _ u . . ., 
T l yield at least n— a— (b— 1) failures. Not all of these failures can occur in the n—a—b 
trials between T n __ a and T 6+1 , respectively. We therefore introduce new events defined for 
0<t<b-l by 

G nt : the n — a— (& — l) st failure among T n _ a , T n „ a „ Xi . . ., T x (in that order) occurs at 2V*. 

Then 

6-1 
Qn-a = 12 G nU SO that 

W.-.)=§ P(Om) for a+b<n<a+c. (49) 

As in the analysis before eq (25), 

W.-«)=S W»i) for a+c<n<2a+c. (50) 

The events F nt and G nt are special cases of the event 
H nkt : the kih failure among T n , T n -i, . . ., T } (in that order) occurs at T, 2 _* +1 _ 
defined for 0<*<&-l and l<Jc<n-b+2. Specifically, 



* 7 nt — rl n -a,n—a—h-\-l, t> F nt — t* n-a,c-~b+l, t- 

Thus it suffices to consider the probabilities P(H nkt ). Now introduce the event 

J n jct'- precisely/:—! failures between T n and T n _ k+2 - t inclusive. 

Then clearly P(H nl0 ) = q n and P(H nkt )=q n _ k +i_ t P(J nkt ) otherwise, so it suffices to consider the 
probabilities P(J nkt ) for (k,t)^(l,Q). 

These probabilities can be computed by recursion on k and t after disposing of the case 
Jc=l by 

P(j nU )=rip n - t iovtyo. 

The initiation of the recursion, for fixed k^>\, is 

P(J„, )=n\ irt _, forifc>l, 
and the "recursive step" is 

P(Jnkt)=Pn- k+ 2-tP(Jn,k,t-l) + qn-^2-tP(Jn,Ic-l,t-l) tOT k>l, t>l. 

7. A Related Problem 

Let K n be the event thai T n is the last of a sequence of a-\-b consecutive trials all yielding 
successes. That is, 

a-\-b-l 

K N = n S n -i. (51) 

For the application which motivated this paper, it is of interest to consider the conditional 
probabilities 

Z h =p(k„IK h [\ £ K) (52) 

defined for n>a-\-b. From the definition of conditional probability, we have 11 

A\,z„-/'(/unA;,n"n' /<:;), 

which because of the inclusion K n ^E n can be rewritten 

X n Z n =P (K n n*nEl). (53) 

Since the preceding material in principle permits the evaluation of X n when a >b, we need only 
investigate the (unconditional) probability on the right-hand side of eq (53). It will be assumed 
that n>a+c+l and &>0. 

It will be shown below that 



which implies that 



K n nK- a =K n n n f\ b s' iy (54) 

i = n— a—c 



n — l n—a—b n—2 

A'„n n E' t =K n n n s;n n e\. 

i=l i = n—a—c i—1 



Recall the definition of X n In eq (2). 
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Since S t ^E i9 the last equation can be rewritten 

7i — l n— a— b 7i—2 7i—a-~c—l 

km n E' t =K n r\ n s;n n em n e\. (55) 

i = l i=n—a—c i=n—a—b+l i=l 

To prove eq (54), first assume its left-hand event occurs. Tentatively suppose Si occurs, 
with n — a—c<i<n—a—b. Since K n occurs, trials T n -\ through T w _ a all yield successes 
(i.e., R n -i holds), and in addition the 6—1 trials T n - a -i through T n _ a _ b+l yield 6—1 successes. 
These (6—1) successes, together with S i} occur between T n _ _i and T w _ a _ c inclusive, and so 
Qn-i-a holds. Thus E n ^i=R n - ) nQ n ^ 1 - u holds, contradicting the assumption that the left- 
hand event in eq (54) holds. So the tentative supposition is untenable, and the left-hand 
event is contained in (i.e., implies) the right-hand one. Conversely, if the right-hand event 
holds then trials T„_ a _ & through T w _ a _ c inclusive constitute (c— 6) + l failures between T n _ a -\ 
and T n ^ a ^ c inclusive, so that Q n _ l _ a , and thus E n _ 1 , must hold. This shows that the right- 
hand event contains the left-hand one, and so eq (54) is proved. 

Next it will be shown that 

n—a—b 

fl S't^JEi lor n-a-b + \<i<n-2. (56) 

i=n—a—c 

This implies that eq (55) can be rewritten 

n—\ n—a—b n—a—c—1 

km n e' ( =km n sm n E' u 

i = l i=n— a — c i = l 

and since the three events on the right-hand side are independent, eq (53) becomes 

/ n—a—b \ /n—a—c—1 \ 

x n z„=i>(K„)P( n s'ApI n e'A 

\i = n—a—c / \ i = l / 

/a+b-1 \ / a+c \ 

=( n Pn-t)( n g[ n AY n - a - c - u (57) 

\ »=0 / \i=a+b / 

thus completing the solution. 

To prove the inclusion (56), suppose its left-hand event occurs, and assume n — a—b-\- 
l<i<n—2. Tentatively assume E t occurs. Then T t through T z _ a +i inclusive must yield 
successes, and so (since i^>n—a—b) must have subscripts greater than the failure-yielding T w _ a _ 6 . 
Thus i— &+l>ft— a— 6, and so n — b<i<n—2. This implies that 

i— a— c+Kjn— a— c<ji— a— b<i— a, 

so that the (c—b) + l failure-yielding trials T n _ a _ & through T n _ a _ c inclusive lie between T^ a 
and 7V a _ c+ i inclusive, thus ensuring that Qi_ a holds. Thus E t holds, contradicting the 
tentative assumption, which is therefore untenable. Therefore (56) is proved. 

Suppose for example that 6 = 1 and c=2. For n>a+c+l=a+3, eq (57) yields 

X n Z n — [ II p n -i ) q n - a -iq n -a-2Y n -a-3. (58) 

\i=0 / 

On the other hand, if 7i>2a+3 then the last equation of section 4, with n replaced by n— 1, 
yields 

<&n—Yn-\ J- n~<ln-a ( II p n -j ) Y n-a-1 ( H Pn-j ) <ln-aQ_n-a-\Y n-a-1 
\j=0 / \y=o / 

+ ( II Pn-j) q n -a-l<Ln-a-2Y n -.a-Z- 

88 



Comparison of this result with eq (58) shows that 

q„_ a ( n Pn-n F n _ a _i — ( II p n -A q n - a qn-a-lYn-a-2=(l—Z n )X n , 

and thus that 

(1— Z w )/Z n ==2»- a (F n _ a _i— q n _ a _ 1 Y n _ a _ 2 )/p n -a<Ln-a-l<Ln-a-2Yn-a-2- 59) 

[f //>3a + 3, however, then eq (27), with n replaced by n—a—l, yields (for a>l) 

o-l /k-l \ 

% n-a-1 Q_n-a-ll n-a-2 = P n-a-l ( [n-a-2~T 2-J Q.n-a-l-k\ H Pn-a-l-j ) J- n-a-k-2 

lc=2 \j=0 J 

+ 1 IT p n - a -]_i ) g n -2a-l<Ln-2a-2Y n-2a-Z- (60) 



On the one hand, eq (60) shows that 

■* n—a Qri—a— 1-* n-a—2cL Pn-a-\^n-a-2^ n-a-3 

so that eq (59) yields 

(1 — Z n )/Z n >2»-a^n-a-l/^n-a2n-a-l. (61) 

This shows, e.g., that if the sequence {p n } is nonincreasing for all 7£>w >2a+3, then Z n <l/2 
for all w>7i +&- A bit less crudely, since the sequence { Y n } is nonincreasing by its definition, 
eq (60) yields 

{a-l /k-\ \ 

( Ln-a-2~\~Z—i Q.n-a-1-k I H Pn-a-l-j ) 

+ ( IT p n _ B _y.iJ 2„-2a-l2»-2a-S f> 

from which a more precise lower bound on (1 — Z„)/Z n can be derived ma eq (59). On the 
other hand, since Y m >q m Y m _i, we have 

' n-a-k-2^1- n-a-3 H ( ]n-a-i> 
I i=3 

/a+2 

F n _2a-3^2 n-a-3 H q n - a - i} 
I i='6 

which when substituted into eq (60) yield an inequality that can be combined with eq (59) to 
obtain a crude upper bound on (1— Z w )/Z w . 

If a=l (as well as 6— -1 and c = 2) then eq (27) yields 

Y n -2 Q71-2-* n-3—Pn -2^n -3<Z« -*¥ »-5j 

so that eq (59) becomes 

(l-Z n )IZ n =q n _ 1 p n . 2 q n ^ i Y n _ 5 /p n _ 1 q n _ 2 Y n _ 4 . (62) 

Since q n ^Y n _ 5 < Y n _ 4 < Y n _ 5 , we have 

qn-iPn-2qn-i/Pn-iqn-2<(l-Z n )/Z n <q tl _ l p n _ 2 /p n _ l q n _ 2 . (63) 

In particular, if {p n } is ultimately nondecreasing then ultimately Z n >l/2. 

(Paper 66B3-75) 
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